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INTRODUCTION 


Specimen  preparation 


The  design  of  structures  in  cold  regions  some¬ 
times  requires  information  on  the  tensile  strength 
of  frozen  soils.  Unfortunately,  very  few  studies 
have  been  done  on  this  subject  to  date.  Haynes  et 
al.  (1975)  investigated  the  effect  of  strain  rate  on 
the  tensile  strength  of  frozen  Fairbanks  silt  with 
high  density  at  -9.4  X.  A  ductile-brittle  transition 
was  reported,  occurring  at  a  strain  rate  of  about 
10'/s. 

The  purpose  of  this  study  is  to  investigate  the  ef¬ 
fect  of  temperature,  strain  rate,  and  density  on  the 
tensile  strength  of  remolded  saturated  frozen  Fair¬ 
banks  silt.  The  test  program  was  as  follows:  1)  the 
samples  with  medium  density  were  tested  under 
various  machine  speeds  from  5.9  x  10'1  to  5.9  x  10’ 
cm/min  at  -2  and  -5X;  2)  the  medium  density 
specimens  were  also  tested  at  various  temperatures 
from  -1  to  -10X  under  a  machine  speed  of  5.9  x 
10"’  cm/min;  and  3)  the  specimens  with  three  dif¬ 
ferent  densities  were  tested  at  the  machine  speed 
of  5.9x10'’  cm/min  and  a  temperature  of  -5°C. 
The  three  desired  densities  are  1.4  (high),  1.2  (me¬ 
dium),  and  1.1  g/cm’  (low),  respectively. 


TESTING 

Material 

The  material  used  in  this  study  was  a  remolded 
silt  taken  from  the  USACRREL  experimental  per¬ 
mafrost  tunnel  at  Fox,  near  Fairbanks,  Alaska. 
The  silt  is  classified  as  ML  in  the  plasticity  chart 
with  approximately  94%  passing  the  #200  mesh 
sieve  and  17%  finer  than  0.01  mm.  Its  main  physi¬ 
cal  properties  are  as  follows: 


Plastic  limit  34.2% 

Liquid  limit  38.4% 

Organic  content  5.49% 

Specific  gravity  2.680 


Molding 

Distilled  water  was  added  to  the  air-dried  Fair¬ 
banks  silt,  making  a  silt-water  mixture  with  an  ini¬ 
tial  molding  water  content  of  12%  by  weight. 
After  being  stored  one  day  to  allow  moisture 
equilibration,  the  moist  silt  was  then  carefully 
compacted  to  the  desired  density  in  a  cast  acrylic 
gang  mold  in  six  layers  as  shown  in  Figure  1 .  The 
mold,  described  in  detail  by  Haynes  (1975),  can 
form  four  25.4-mm-diameter  by  82.55-mm-long 
dumbbell  specimens  at  a  time. 

Saturating 

After  compaction,  the  samples  were  saturated 
with  de-aerated,  distilled  water  under  a  vacuum  of 
73  mm  Hg.  The  saturation  system  is  shown  in  Fig¬ 
ure  2. 

Freezing 

When  the  samples  were  saturated,  they  were 
placed  in  a  freezing  cabinet  and  quickly  frozen 
from  the  top  down  in  an  open  system  at  -30X.  Fig¬ 
ure  3  shows  the  freezing  cabinet  with  a  specimen- 
charged  mold  ready  to  be  set  up  for  freezing.  To 
obtain  rapid  frost  penetration  and  thus  minimize 
ice  lensing  during  freezing,  the  specimens  were  al¬ 
lowed  to  equilibrate  at  OX  in  the  cabinet  over¬ 
night  before  freezing.  Some  minor  ice  lenses  paral¬ 
lel  to  the  frost  penetration  are  still  visible  in  the 
specimens  with  low  density,  even  though  the  freez¬ 
ing  process  is  carefully  controlled.  After  freezing, 
the  mold  was  taken  from  the  cabinet  and  carefully 
disassembled  in  a  coldroom,  as  shown  in  Figure  4. 

Dimension  of  specimen 

One-half  of  the  specimen  geometry  is  illustrated 
in  Figure  5.  The  neck -down  portion  of  the  dumb¬ 
bell-shaped  specimen  is  50.8  mm  long  and  25.4 
mm  in  diameter. 
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Testing  procedure  and  apparatus 

The  tension  tests  were  conducted  on  a  Model 
16000  MTS  closed-loop  electrohydraulic  testing 
machine  operated  in  “stroke”  control.  This  mode 
of  control  is  equivalent  to  a  constant  cross-head 
speed  test.  A  Ransco  refrigeration  unit  with  a 
dual-chambered  temperature-controlled  cabinet 
was  installed  in  conjunction  with  the  MTS  ma¬ 
chine.  The  two  chambers  (an  air-temperature  con¬ 
ditioning  and  a  test  chamber)  are  connected 
through  large  insulated  air  transfer  hoses.  Figure  6 
shows  the  MTS  machine  and  the  installed  testing 
cabinet.  Before  testing,  each  specimen  was  tem¬ 
pered  in  the  cabinet  at  the  test  temperature  for  at 
least  24  hours.  The  environmental  temperature 
around  the  specimen  was  measured  by  a  digital 
thermometer  with  a  sensitivity  of  0.01  °C.  Obser¬ 
vation  has  shown  that  its  fluctuation  is  less  than 
±0.05‘t  during  tests. 

The  axial  tensile  load  was  measured  by  a  multi¬ 
range  load  cell  accessory  to  the  MTS  machine, 
Model  661.21A-03.  The  elongation  of  the  neck 
section  of  a  specimen  was  measured  by  two  exten- 
someters,  which  were  held  to  the  surface  of  the 
specimen  by  rubber  bands,  as  shown  in  Figure  7. 
The  gauge  length  of  the  extensometers  is  50.8  mm 
with  a  10.2  mm  range.  Two  edges  of  each  exten- 
someter  were  carefully  set  up  on  the  surface  at 
both  ends  of  the  neck  section  of  a  specimen.  The 
average  values  of  the  elongation  measured  by  the 
two  extensometers  were  then  used  to  calculate  the 
strain  of  the  specimen.  For  some  of  the  tests,  the 
elongation  of  the  specimens  between  the  two  end 
caps  was  also  measured  by  two  LVDTs  (linear  var¬ 
iable  differential  transformers),  Model  200-DCD, 
which  were  mounted  on  the  end  caps,  as  shown  in 
Figure  8. 

To  protect  specimens  from  sublimation  during 
testing,  the  surface  of  specimens  tested  at  machine 
speeds  less  than  5.9  x  10”  cm/min  was  coated  with 
a  thin  layer  of  silicone  grease. 

During  a  test,  the  tensile  load  and  the  elonga¬ 
tion  of  the  specimen  were  recorded  by  a  strip-chart 
recorder  and  a  Kaye  Model  8400  data  logger  (for 
slow  tests)  or  a  Biomation  Model  1015  transient 
recorder  together  with  a  Type  203  oscilloscope 
(for  fast  tests).  These  two  instruments  are  shown 
in  Figure  9.  The  environmental  temperature  was 
recorded  by  a  strip-chart  recorder. 

After  testing,  the  neck  section  of  each  specimen 
was  removed  for  water  content  and  density  meas¬ 
urement. 


RESULTS 

The  results  of  the  constant  strain-rate  tension 
test  are  summarized  in  Tables  1,  2,  and  3.  The 
peak  (or  ultimate)  strength  om  refers  to  the  maxi¬ 
mum  stresses  on  the  true  stress  vs  true  strain* 
curves.  Correspondingly,  the  failure  strain  tf  and 
the  time  to  failure  fm  refer  to  the  strain  and  time  at 
the  peak  strength,  respectively.  The  initial  tangent 
modulus  E[  and  the  50%  strength  modulus  £,  are 
determined  graphically.  The  average  strain  rate  is 
computed  by  c  =  e<-/fm.  If  there  was  no  data  for  tf 
and  tm,  the  average  strain  rate  was  then  taken  as  a 
nominal  value  corresponding  to  the  applied  ma¬ 
chine  speed. 

The  physical  properties  of  the  specimens  are 
shown  in  Appendix  A.  Note  that  water  content  w 
and  dry  density  7d  are  taken  from  the  specimens 
after  testing.  The  volume  ice  content  /„  is  calculat¬ 
ed  by 

(w-wubd 


where  Gj  is  the  specific  gravity  of  ice  taken  as 
0.917,  yw  is  the  unit  weight  of  water  taken  as  1 
g/cm\  and  wu  is  the  unfrozen  water  content  of  the 
specimens.  The  values  of  wu  of  Fairbanks  silt  with 
various  total  water  contents  are  shown  in  Appen¬ 
dix  B.  The  saturation  degree  S  is  computed  by 

100[w-(l  -Gi)wJG 

S  =  - ^ - *  -  (%)  (2) 

G,e 

where  G  is  the  specific  gravity  of  the  soil  particles, 
equal  to  2.68  for  the  soil  tested,  and  e  is  the  void 
ratio  of  the  specimens. 


DISCUSSION 

Stress-strain  curves  and  failure  mode 

Figure  10  presents  a  set  of  typical  stress-strain 
curves  of  the  samples  with  medium  density  for 


•In  this  study,  the  true  strain  is  calculated  by  (  =  Inf  I  +  tc), 
where  <c  =  1L/L,  is  the  conventional  strain,  in  which  A L  and 
L,  are  the  axial  elongation  and  the  initial  gauge  length  of  speci¬ 
mens.  The  true  stress  is  calculated  by  a  =  P(1  +  (J/A,,  where 
P  is  the  total  axial  tensile  force  and  A,  is  the  initial  cross- 
sectional  area  of  the  specimens. 
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various  strain  rates  at  -S'C.  It  is  seen  that  strain 
rate  has  a  strong  effect  on  the  stress-strain  behav¬ 
ior  of  the  frozen  silt.  The  brittle  behavior  of  the 
frozen  soil  increases  with  increasing  strain  rate.  At 
strain  rates  greater  than  about  10'Vs,  the  soil  exhi¬ 
bits  a  brittle  type  fracture,  which  is  characterized 
by  a  smaller  failure  strain  (less  than  0.5%)  and  a 
nearly  straight  a-t  curve.  For  this  type  of  failure, 
the  elasticity  of  the  material  essentially  governs  its 


mechanical  behavior.  For  i  <  10'2/s  the  sample 
exhibits  a  ductile  or  a  moderately  brittle  type  fail¬ 
ure,  characterized  by  a  significant  plastic  flow  fol¬ 
lowing  a  small  elastic  deformation. 

It  should  be  noted  that  for  relatively  low  strain 
rates  (say,  <  10' Vs),  the  a-t  curves  rapidly  drop 
down  from  their  peak  at  a  small  strain  of  less  than 
0.5%.  This  type  of  a-t  behavior  is  completely  dif¬ 
ferent  from  that  observed  in  compression,  where 


Table  1.  Tension  test  results  of  frozen  Fairbanks  silt  at  -5°C. 


Specimen 

no. 

Machine 

speed 

V 

(cm/min) 

Peak 

strength 

(kg/cm1) 

Failure 

strain 

•J 

<*) 

Time  to 
failure 

fs) 

A verage 
strain  rate 

t 

Us) 

Initial  tangent 
modulus  E, 

(  x  10') 

(kg  /cm1) 

50%  peak  strength 
modulus  E, 
(x)O') 
(kg/cm1) 

12-4 

5.9  x  10' 

58.9 

0.240 

Medium  density 

7.00x10"  3.43x10" 

12-3 

5.9  x  10' 

58.2 

0.390 

8.00  x  10" 

4.88  x  10" 

— 

— 

11-4 

5.9x  10' 

72.5 

— 

8.00x10" 

4.88  x  10" 

— 

— 

12-2 

5.9  x  10' 

63.6 

— 

8.00x10" 

4.88x10" 

— 

— 

11-3 

5.9x10' 

73.4 

— 

8.00x10" 

4.88  x  10" 

— 

— 

12-1 

5.9x  10" 

65.5 

0.770 

7.80  x  10" 

9.87x10" 

— 

— 

5-2 

5.9x10" 

81.1 

0.650 

7.00x10" 

9.30x10" 

— 

— 

6-1 

5.9  x  10" 

75.2 

0.800 

1.00x10* 

8.00x10" 

— 

— 

4-2 

5.9x10" 

44.4 

1.215 

1.20x  10' 

1.01x10" 

— 

— 

4-3 

5.9x10" 

48.7 

1.500 

1.50x10' 

1.00x10" 

— 

— 

17-2 

5.9x10-' 

47.1 

1.824 

1.90x10’ 

9.60x10" 

— 

— 

17-4 

5.9x10-' 

51.6 

1.482 

1.40x10' 

1.06x10" 

— 

— 

3-4 

5.9x10-' 

31.6 

1.049 

l.lOxlO1 

1.00x10" 

10.0 

— 

4-1 

5.9x10’ 

32.0 

1.128 

l.lOx  1C 

1.03x10" 

11.0 

1.35 

4-4 

5.9x10-’ 

22.4 

1.224 

1.24X  10’ 

1.07x10" 

5.0 

0.93 

5-3 

5.9x10" 

21.0 

0.574 

7.24  x  10’ 

7.93  x  10" 

13.0 

1.00 

5-4 

5.9x10" 

22.6 

0.936 

9.04X1C 

1.04x10" 

12.0 

1.30 

6-2 

5.9x10“" 

18.0 

0402 

4.86x10’ 

8.27  x  10" 

14.8 

1.80 

6-3 

5.9x10" 

19.2 

0.477 

6.30x10’ 

7.57x10" 

11.0 

2.00 

11-1 

5.9x10-’ 

16.7 

0.334 

4.44  x  1C 

7.52x10" 

9.5 

1.58 

11-2 

5.9x10-’ 

14.8 

0.295 

4.68  x  1C 

6.30  x  10" 

9.0 

1.30 

15-1 

5.9x10  ' 

44.3 

2.608 

High  density 

2.40x10' 

1.09x10" 

15-2 

5.9x10-' 

43.7 

2.585 

2.40x10' 

1.08x10" 

— 

— 

9-2 

5.9x10-' 

27.8 

1.665 

1.52x10’ 

1.10x10" 

12.0 

0.80 

9-1 

5.9x10-' 

29.8 

1.946 

1.40x10' 

1.39x10" 

12.0 

0.89 

15-3 

5.9x10-’ 

20.3 

1.282 

1.46x10’ 

8.78x10" 

12.0 

1.00 

9-3 

5.9x10-' 

25.3 

1.660 

1.38x10’ 

1.20x10" 

9.5 

1.10 

9-4 

5.9x10* 

16.4 

0.387 

5.64x10’ 

6.86x10" 

9.0 

1.30 

15-4 

5.9x10- 

15.3 

0.698 

7.32x10' 

9.54x10" 

6.8 

0.81 

18-2 

5.9x10* 

41.2 

0.250 

Low  density 

4.00x10" 

6.25x10" 

18-3 

5.9x10* 

45.7 

— 

8.00x10" 

6.25x10" 

— 

— 

18-1 

5.9x10 " 

51.9 

0.480 

1. 00x10' 

4.80x10" 

— 

— 

16-1 

5.9x10" 

46.5 

0.389 

7.00  x  1C 

5.56x10" 

— 

— 

10-1 

5.9x10" 

39.9 

0.998 

2.02x10’ 

4.95  x  10" 

13.7 

1.54 

10-2 

5.9x10 -’ 

37.1 

0.874 

1.88x10’ 

4.65  x  10" 

11.0 

1.50 

10-3 

5.9x10 

26.7 

0.586 

1.35x10’ 

4.33  x  10" 

19.0 

1.40 

16-2 

5.9x10" 

25.4 

0.632 

1 .02  x  10’ 

6.22x10" 

18.3 

1.36 

10-4 

5.9x10" 

20.4 

0.548 

1 .06  x  1C 

5.19x10" 

16.8 

1.60 

16-3 

5.9x10" 

17.5 

0.449 

8.00x10’ 

5.61  x  10" 

13.0 

1.40 
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Table  2.  Tension  test  results  of  frozen  medium-density  Fairbanks  silt  at  -2°C. 


Specimen 

no. 

Machine 

speed 

V 

( cm/min ) 

Peak 

strength 

( kg/cm ’) 

Failure 

strain 

7 

(%) 

Time  to 
failure 
t/n 
(s) 

A verage 
strain  rate 

t 

Us) 

Initial  tangent 
modulus  E, 
(*10‘) 
(kg/cm1) 

50%  peak  strength 
modulus  E, 
(xlO-J 
(kg  /cm1) 

19-3 

5.9x10' 

40.7 

9.00x10" 

19-4 

5.9x10' 

38.8 

— 

— 

9.00x10" 

— 

19-1 

5.9x10“ 

48.4 

0.540 

6.00x10" 

9.00x10" 

— 

— 

19-2 

5.9x10“ 

52.4 

— 

9.00x10" 

— 

— 

13-2 

5.9x10-' 

30.4 

1.987 

2.00  - 10 

1.00x10" 

— 

— 

13-3 

5.9x10" 

26.1 

1.751 

1.60x10' 

1.60x10" 

— 

— 

13-4 

5.9x10-' 

26.6 

1.550 

1.60x10' 

9.69x10" 

— 

— 

8-3 

5.9  x  10" 

20.9 

1.465 

1.85x10' 

7.92x10" 

6.0 

0.84 

8-4 

5.9X10" 

23.6 

1.301 

1.32x10’ 

9.86x10" 

7.0 

0.80 

2-1 

5.9x10" 

19.0 

1.290 

1.68x10' 

7.69x10" 

6.5 

— 

1-3 

5.9x10-' 

14.6 

0.858 

3.60x10' 

2.38x10" 

6.9 

0.70 

14-3 

5.9x10" 

7.8 

0.900 

8.64x10* 

1.04x10" 

4.6 

0.35 

14-4 

5.9x10" 

7.7 

0.881 

8.64x10* 

1.02x10" 

5.0 

0.30 

Table  3.  Tension  test  results  of  frozen  medium-density  Fairbanks  silt  at  a  machine  speed  of 
0.059  cm/min. 


Specimen 

no. 

Temperature 

e 

(°C) 

Peak 

strength 

(kg /cm’) 

Failure 

strain 

7 

(%) 

Time  to 
failure 

(S) 

A verage 
strain  rate 

t 

Us) 

Initial  tangent 
modulus  E, 
(xlO') 
(kg/cm1) 

50%  peak  strength 
modulus  E, 
(xIO'/ 
(kg/cm’) 

14-1 

-  1.0 

16.8 

2.194 

2.20x10’ 

9.97x10" 

6.0 

0.55 

14-2 

-  1.0 

16.1 

2.665 

2.72x10’ 

9.80x10" 

5.6 

0.50 

8-3 

-  2.0 

20.9 

1.465 

1.85x10' 

7.92x10" 

6.0 

0.84 

8-4 

-  2.0 

23.6 

1.301 

1.32x10’ 

9.86x10" 

7.0 

0.80 

2-1 

-  2.0 

19.0 

1.290 

1.68x10' 

7.69x10" 

6.5 

0.80 

8-1 

-  3.0 

25.2 

1.153 

1.20x10’ 

1.20x10" 

8.5 

1.00 

8-2 

-  3.0 

25.2 

1.533 

1.48x10’ 

1.04x  10" 

8.5 

1.00 

3-4 

-  5.0 

31.6 

1.049 

1.10x10’ 

1.00x10" 

10.0 

— 

4-1 

-  5.0 

32.0 

1.128 

1.10x10’ 

1.03x10" 

11.0 

1.35 

7-3 

-  7.0 

41.1 

1.257 

1.50x10’ 

8.38  x  10" 

12.8 

1.80 

7-4 

-  7.0 

41.5 

1.028 

1.40x10’ 

7.34  x  10" 

14.0 

1.80 

6-4 

-10.0 

58.6 

0.837 

1.30x10’ 

6.44  x  10" 

16.0 

2.40 

7-1 

-10.0 

50.9 

1.468 

1.86x10' 

7.90x10" 

17.6 

2.35 

7-2 

-10.0 

56.8 

1.209 

1.52x10’ 

7.95x10" 

16.8 

2.20 

the  materials  sustain  a  significant  stress  level  after 
initial  yielding.  Because  of  this,  great  care  must  be 
taken  in  engineering  design  where  the  long-term 
tensile  strength  of  frozen  soils  needs  to  be  consid¬ 
ered. 

Figure  1 1  shows  a  set  of  o-t  curves  of  the  medium- 
density  samples  for  various  temperatures  at  a  strain 
rate  of  10Vs.  Obviously,  the  brittleness  of  the  sam¬ 
ple  increases  with  decreasing  temperature.  However, 
the  failure  mode  is  less  dependent  on  temperature 
than  on  strain  rate.  All  of  the  o-t  curves  for  temper¬ 
atures  varying  from  -1  to  -1°C  (Fig.  11)  reflect  a 


similar  level  of  plasticity.  Even  at  a  lower  tempera¬ 
ture  of  - 1 0 °C,  no  brittle  fracture  was  observed. 

Figure  12  presents  some  typical  o-t  curves  of  the 
specimens  with  different  dry  densities  at  a  strain  rate 
of  about  IO'Vs  and  a  temperature  of  -S°C.  The  o-t 
curve  appears  to  be  somewhat  sharper  for  the  speci¬ 
men  with  lower  density  than  for  higher  density,  indi¬ 
cating  a  slight  loss  of  ductility  with  decreasing  den¬ 
sity. 

Two  specimens  after  testing  are  shown  in  Figure 
13.  It  seems  that  as  brittleness  increases,  the  fracture 
plane  of  the  specimens  becomes  flatter. 
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Peak  tensile  strength,  am 

Strain  rate  dependence  of 
the  peak  strength 

The  peak  tensile  strength  am  of  the  frozen  silt  with 
medium  density  as  a  function  of  i  at  -2  and  -5^  is 
plotted  in  Figure  14.  The  ductile-brittle  (D-B)  tran¬ 
sition  occurs  at  almost  the  same  strain  rate  of  about 
10'Vs  at  both  -2  and  -5^,  indicating  that  the  criti¬ 
cal  strain  rate  of  the  D-B  transition  is  not  sensitive 
to  temperature  in  this  temperature  range  for  the  me¬ 
dium-density  Fairbanks  silt.  It  was  found,  however, 
that  this  critical  strain  rate  varies  with  dry  density. 

Figure  15  presents  a  log-log  plot  of  am  vs  e  for  var¬ 
ious  dry  densities  at  -5^.  It  shows  that,  for  lower 
density,  the  critical  strain  rate  decreases  to  about  5  x 
10'Vs.  Based  on  this  observation,  one  may  infer  that 
the  critical  rate  of  D-B  transition  for  ice  should  be 
lower  than  the  above  value.  According  to  Mellor 
(1972),  this  is  true.  He  reported  that  the  ductile-brit¬ 
tle  transition  occurs  at  a  strain  rate  of  about  10~4/s 
for  polycrystalline  ice  at  -7^.  Unfortunately,  we 
did  not  run  the  high  strain  rate  tests  on  the  high-den¬ 
sity  samples,  so  the  D-B  transition  rate  for  high  den¬ 
sity  could  not  be  determined  in  this  study.  According 
to  Haynes  (1975),  it  is  about  10'Vs  at  -9.4C'C. 

Note  that  Figure  14  shows  another  break  point  at 
a  strain  rate  of  10'Vs  on  the  logt  -  log<rm  curve  for 
-S'C.  Reviewing  Figure  10,  one  can  see  that  a  signif¬ 
icant  drop  in  failure  strain  (from  about  1%  to  0.4%) 
also  occurs  at  this  strain  rate.  This  evidence  may  in¬ 
dicate  a  change  in  deformational  mechanisms  at  this 
critical  strain  rate. 

From  Figures  14  and  15,  the  peak  strength  om  as  a 
function  of  strain  rate  can  be  described  by  a  simple 
power-law  equation  within  the  range  of  strain  rate 
applied: 

°m  =  *(f/f.)"  (3) 

where  f,  is  a  reference  strain  rate  taken  as  1.0/s,  n  is 
a  dimensionless  parameter,  and  k  has  the  dimension 
of  stress.  Both  k  and  n  vary  with  temperature,  densi¬ 
ty,  and  the  range  of  strain  rates.  The  values  of  k  and 
n  determined  from  the  test  data  by  linear  regression 
analysis  are  shown  in  Table  4. 

It  is  interesting  to  compare  the  strain  rate  depen¬ 
dence  of  tensile  strength  with  that  of  compressive 
strength.  The  comparison  for  -5°C  is  shown  in  Fig¬ 
ure  16,  in  which  the  compression  data  are  from  the 
authors’  recent  report  (Zhu  and  Carbee  1984).  It  is 
seen  from  this  figure  that  for  i less  than  10'!/s,  the 
peak  tensile  strength  is  close  to  the  peak  compressive 
strength.  However,  for  c  greater  than  10'Vs,  the 


Table  4.  Values  of  k  and  n  in  equation  3. 


f) 

ro 

Id 

fg/ cm') 

t 

f/s) 

k 

(kg/cm‘) 

n 

1.08-1.12 

5x  10  *  -  5  x  10  ’ 

148.6 

0.142 

-5 

1.20-1.26 

1  x  10  '  -  1  x  10’’ 

1  x  10  ’  -  6x  10  • 

143.4 

48.6 

0.151 

0.068 

1.36-1.41 

1  x  10  ’  -  7x10-’ 

105.4 

0.134 

-2 

1.20-1.26 

1  x  10-'  -  1  x  10  * 

103.2 

0.185 

former  is  much  less  than  the  latter.  Haynes  (1975) 
observed  a  similar  bifurcation  at  a  strain  rate  of 
about  10'Vs  for  Fairbanks  silt  with  high  density  at 
-9. 4^0.  Hawkes  and  Mellor  (1972)  reported  a  bifur¬ 
cation  of  the  two  strengths  at  about  10'Vs  for  poly¬ 
crystalline  ice.  Further  micromechanistic  studies  are 
needed  to  discover  the  physical  reasons  for  this  bi¬ 
furcation. 

Time  dependence  of  the  peak  strength 

The  peak  tensile  strength  om  as  a  function  of 
time  to  failure  tm  for  various  dry  densities  and  tem¬ 
peratures  is  shown  in  Figures  17  and  18,  respec¬ 
tively.  It  is  obvious  that  the  peak  strength  as  a 
function  of  time  to  failure  tm  has  an  upper  limit, 
i.e. ,  the  limiting  instantaneous  strength  a„. 

In  a  ductile  failure  regime,  am  rapidly  decreases 
with  increasing  ?m  until  tm  reaches  about  200  s, 
and  then  am  decreases  very  slowly  with  increasing 
tm.  In  a  brittle  fracture  regime,  am  decreases 
slightly  with  decreasing  /m.  The  average  values  of 
o„  obtained  from  the  present  study  are  shown  in 
Table  5.  The  time  to  failure  corresponding  to  a0  is 
relatively  insensitive  to  temperature  for  the  range 
tested,  but  varies  with  density.  It  has  a  value  of 
about  1  second  for  the  silt  with  medium  density 
and  about  10  seconds  for  low  density. 

Replotting  Figure  18  in  the  coordinates  of  l/om 
vs  log/m,  as  shown  in  Figure  19,  one  can  see  that 
the  tensile  strength  of  the  frozen  silt  as  a  function 
of  time  to  failure  for  the  ductile  failure  regime  can 


Table  5.  Average  values  of  o,  for 
frozen  Fairbanks  silt. 


o0  (kg/cm') 


V 

Medium 

Low 

rcj 

density 

density 

-2 

50.4 

_ 

-5 

73.9 

49.2 

be  well  described  by  Vialov’s  (1959)  strength  equa¬ 
tion: 


Table  7.  Values  of  A  and  m  in  equation 
5  for  (  =  10'Vs  and  =  1.20-1.26 
g/cm\ 


_ 0 _ 

am  log(/m/B) 


(4) 


0  A 

(°C) _ (kg  /cm1)  m  7  N 


where  0  and  B  are  empirical  parameters.  0  has  the 
dimension  of  stress,  and  B  has  the  dimension  of 
time.  For  a  particular  type  of  soil,  both  0  and  B 
depend  upon  temperature  and  density.  The  values 
of  0  and  B  obtained  from  the  test  data  are  shown 
in  Table  6.  Note  that  the  strength-time  data  for 
brittle  fracture  could  not  be  used  to  determine  the 
values  of  0  and  B  in  eq  4. 

Temperature  dependence  of 
the  peak  strength 

Figure  20  presents  a  log-log  plot  of  am  vs  6/60 
for  the  samples  with  medium  density  at  an  average 
strain  rate  of  about  10'Vs.  It  shows  a  break  point 
at  the  temperature  of  -5°C,  and  the  peak  strength 
increases  more  rapidly  with  decreasing  tempera¬ 
ture  when  0  <  -5^.  From  this  figure,  am  as  a 
function  of  d  can  be  described  by  a  power-law 
equation  within  a  certain  range  of  temperature: 

am  =  A  (0/0.r  (5) 

where  0  is  the  negative  temperature  in  °C,  60  is  a 
reference  temperature  taken  as  -1  °C,  and  A  and  m 
are  empirical  parameters.  From  these  test  results, 
they  depend  upon  temperature.  The  values  of  A 
and  m,  together  with  the  correlation  coefficient  7 
and  data  points  N  obtained  from  the  test  results, 
are  shown  in  Table  7. 

Figure  21  presents  a  log-log  curve  of  unfrozen 
water  content  against  temperature  for  the  silt  with 
a  water  content  of  40.5%,  which  is  plotted  using 
the  data  shown  in  Appendix  B.  It  is  interesting  to 
note  that  the  log-log  curve  of  wu  vs  0/6o  also 
breaks  at  -S^C  (Fig.  21). 


Table  6.  Values  of  0  and  B  in 
equation  4. 

e  Id  fl  B 


ro 

(g/cm') 

(kg/cm') 

(s) 

1.08-1.12 

66.7 

3.4 

-5 

1.20-1.26 

80.1 

0.3 

1.36-1.41 

65.9 

1.1 

-2 

1.20-1.26 

21.4 

14.6 

-1  to  -5  16.5  0.404  0.9892  8 

-5  10  -1  8.7  0.801  0.9855  7 


Effect  of  density  on  the  peak  strength 

The  peak  strength  am  as  a  function  of  dry  den¬ 
sity  7d  for  various  machine  speeds  at  -5°C  is  plot¬ 
ted  in  Figure  22.  It  is  seen  that  for  the  saturated 
frozen  silt,  am  increases  with  decreasing  7^  in  the 
ductile  regime.  Especially  for  relatively  high  strain 
rates  and  low  density  range,  am  is  more  sensitive  to 
7<j.  As  is  expected,  this  relationship  can  also  be 
seen  from  a  plot  of  am  vs  as  shown  in  Figure  23, 
in  which  the  volume  ice  content  iv  is  calculated  by 
eq  1 .  A  similar  relationship  was  observed  between 
an  initial  compressive  yield  strength  and  dry  densi¬ 
ty  or  between  that  and  volume  ice  content  (Sayles 
and  Car  bee  1981,  Zhu  and  Carbee  1984),  indicat¬ 
ing  that  the  tensile  strength,  like  the  initial  yield 
strength  in  compression,  is  essentially  governed  by 
the  fracturing  of  the  ice  matrix  in  frozen  soil. 

Failure  tensile  strain,  ef 

As  defined  before,  the  failure  tensile  strain  ef  re¬ 
fers  to  the  strain  at  the  peak  tensile  strength  am. 
Figure  24  shows  failure  strain  ef  as  a  function  of  e 
for  the  specimens  with  various  densities  at  -5°C. 
Explicitly,  ef  increases  with  increasing  e  in  the  duc¬ 
tile  failure  regime,  while  it  decreases  with  increas¬ 
ing  e  in  the  brittle  regime.  The  maximum  failure 
strain  occurs  at  a  strain  rate  of  about  10'Vs  for 
medium  density  and  5  x  10'Vs  for  low  density, 
which  is  one  order  less  than  that  corresponding  to 
the  limiting  instantaneous  strength  a0. 

Figure  24  also  shows  that  for  relatively  high 
strain  rates,  failure  strain  increases  considerably 
with  increasing  dry  density.  For  example,  at  a 
strain  rate  of  10' Vs,  e f  increases  from  about  0.4% 
to  2.6%  as  7d  increases  from  1.1  to  1.4  g/cm’  at 
-5°C.  However,  there  is  not  a  significant  differ¬ 
ence  between  the  failure  strains  of  the  specimens 
with  different  densities  at  lower  strain  rates. 

The  failure  strain  as  a  function  of  temperature 
for  a  given  machine  speed  and  density  is  shown  in 
Figure  25.  It  is  clear  that  ef  rapidly  increases  with 
increasing  temperature  for  temperatures  higher 
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than  -2°C.  This  may  be  explained  by  the  fact  that 
as  temperature  rises  above  the  significant  phase 
transition  temperature  of  about  -2°C  for  the  silt, 
the  rapid  increase  in  the  unfrozen  water  content 
with  temperature  greatly  increases  the  plasticity  of 
the  soil. 

The  failure  strains  from  both  the  tension  and 
compression  tests  on  the  medium-density  speci¬ 
mens  are  compared  in  Figure  26.  shows  a  similar 
variation  with  i  for  both  tension  and  compression. 
However,  the  magnitude  of  t f  in  the  tension  tests  is 
about  one  order  less  than  in  the  compression  tests. 

Initial  tangent  modulus  and 
50%  peak  strength  modulus 

The  initial  tangent  modulus  E\  and  modulus  £, 
at  50%  strength  of  the  medium-density  samples  as 
a  function  of  temperature  are  plotted  on  log¬ 
arithm  coordinates  in  Figures  27  and  28,  respec¬ 
tively.  The  curves  in  the  two  figures  break  at 
-5°C;  they  can  be  described  by  a  power-law  equa¬ 
tion  within  the  range  of  temperature  tested: 


£;,£,  =  a(0/6o)b 


(6) 


For  the  sake  of  comparison,  the  E\  data  of  frozen 
Fairbanks  silt  reported  by  various  researchers  are 
cited  as  follows:  Kaplar  (1963)  reported  Young’s 
modulus  values  from  2.25x10’  to  2.39  x 
10’  kg/cm2  within  a  range  of  temperature  from  0 
to  -10°C  by  a  dynamic  method.  Offensend  (1966) 
found  an  average  value  of  2  x  10’  kg/cm2  from  his 
tension  tests  at  both  -3.9  and  -9.4 °C;  and  Haynes 
(1975)  reported  an  average  value  of  about  7  x  10* 
kg/cm2  from  his  tension  tests  on  dense  samples  at 
-9.4°C. 

Discussion  of  test  techniques 

End  caps 

The  experiment  shows  that  most  of  the  samples 
tested  under  intermediate  and  slow  speeds  failed 
within  gauge  length.  However,  the  specimen  tested 
under  fast  speeds  (say,  >  2xI0"2/s)  frequently 
broke  near  the  end  cap,  as  shown  in  Figure  32. 
This  is  probably  due  to  the  stress  concentrating  in 
the  specimen  near  the  end  caps  during  loading. 
Further  studies  need  to  be  done  to  improve  the  end 
cap  configuration. 


where  parameters  a  and  b  depend  upon  the  tem¬ 
perature  range  for  a  given  soil.  The  values  of  a  and 
b  obtained  from  the  test  data  are  given  in  Table  8. 

The  initial  tangent  modulus  and  the  50% 
strength  modulus  as  a  function  of  i  for  various 
densities  at  -5°C  are  shown  in  Figures  29  and  30, 
respectively.  It  is  seen  that  both  E\  and  £,  are  not 
dependent  upon  «,  but  tend  to  increase  with  de¬ 
creasing  dry  density. 

Figure  31  presents  a  plot  of  E,  against  £,  for 
various  test  conditions.  There  seems  to  be  a  defi¬ 
nite  relationship  between  E\  and  £,,  even  though 
the  data  are  scattered;  it  can  be  expressed  as 


£j  =  6.05£,  +  3.26x10* 


(7) 


where  E\  and  £,  are  in  kg/cm2.  This  relationship  is 
helpful  for  determining  £j  because  £,  is  easier  to 
determine  from  the  test  data. 


Elongation  measurement 

Measuring  the  elongation  of  the  neck  section  of 
a  specimen  in  tension  tests  is  a  difficult  problem. 
Considering  the  geometric  configuration,  Haynes 
(1975)  presented  a  reduction  procedure  for  calcu¬ 
lating  the  strain  in  the  neck  section  of  a  dumbbell 
specimen,  based  on  the  measured  deformation 
between  the  two  end  caps.  Following  his  proce¬ 
dure,  we  deduced 


AL  =  0.665AL, 


(8) 


where  A Lx  is  the  observed  deformation  of  the  en¬ 
tire  specimen  between  the  two  end  caps  and  A L  is 
the  corrected  deformation  in  the  neck  section. 

In  the  present  study,  the  deformation  in  the  neck 
section  of  a  specimen  was  directly  measured  by  us¬ 
ing  two  extensometers.  As  discussed  earlier,  for 
the  sake  of  comparison,  the  deformation  of  the 


Table  8.  Values  of  a  and  b  in  equation  6. 

E,  E, 


0 

ro 

a 

Cx/OV 

(kg/cm’) 

b 

> 

l V 

a 

(x!0‘) 

(kg/cm’J 

b 

1 

N  1 

& 

I' 

01  to  05 

5.50 

0.381 

0.9460 

9 

0.531 

0.586 

0.9946 

8  I 

?v 

05  to  010 

3.56 

0.676 

0.9758 

7 

0.402 

0.762 

0.9897 
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entire  specimen  between  the  two  end  caps  was  also 
measured  for  some  of  the  tests,  and  eq  8  was  used 
to  compute  the  deformation  in  the  neck  section  of 
the  specimen.  The  observed  and  the  calculated  val¬ 
ues  of  the  strains  at  failure  are  compared  in  Figure 
33.  It  shows  that  the  strain  up  to  failure  calculated 
by  using  Haynes’  reduction  procedure  is  slightly 
(about  6%)  less  than  the  measured  values. 

Prevention  of  sublimation 
The  specimen  data  in  Table  B1  show  that  the 
saturation  degree  of  the  moist  specimens  coated 
with  silicone,  determined  after  the  tests,  is  close  to 
100%.  That  is  to  say,  coating  a  specimen  with  a 
thin  layer  of  silicone  does  protect  the  specimen 
from  sublimating  for  a  certain  period  of  elapsed 
time.  To  determine  whether  the  silicone  has  an  ef¬ 
fect  on  the  strength  of  the  specimen,  the  following 
tests  were  carried  out.  Two  replicate  specimens 
were  carefully  prepared.  One  was  enclosed  in  a 
plastic  bag  with  the  plastic  mold  inserts  (used  to 
form  the  necked  portion  of  the  specimen)  still  in 
place.  The  other  was  coated  with  silicone.  After 
the  two  samples  were  tempered  in  an  environmen¬ 
tal  chamber  for  48  hours,  they  were  tested  at  the 
same  conditions.  As  shown  in  Figure  34,  very  sim¬ 
ilar  a-e  curves  were  observed  for  these  two  sam¬ 
ples,  indicating  that  the  silicone  layer  has  no  sig¬ 
nificant  effect  on  the  results. 


CONCLUSIONS 

This  investigation  leads  to  the  following  conclu¬ 
sions. 

1 .  The  transition  from  ductile  failure  to  brittle 
failure  was  observed  occurring  at  a  certain  strain 
rate.  This  critical  strain  rate  was  found  not  to  be 
sensitive  to  temperature  for  temperatures  down  to 
-5°C,  but  it  varies  with  density.  It  is  about  10'!/s 
for  medium-density  and  5  x  10"4/s  for  low -density 
silt  at  -5°C.  For  brittle  fracture,  the  peak  tensile 
strength  slightly  decreases  with  increasing  strain 
rate,  while  for  ductile  failure,  it  significantly  de¬ 
creases  with  decreasing  strain  rate. 

2.  The  peak  tensile  strength  increases  with  de¬ 
creasing  temperature,  and  it  increases  more  rapid¬ 
ly  when  the  temperature  is  lower  than  -5°C. 

3.  The  peak  tensile  strength  increases  with  de¬ 
creasing  density,  especially  for  lower  densities. 


4.  The  failure  strain  varies  with  strain  rate  and 
density,  but  does  not  depend  upon  temperature 
when  temperature  is  lower  than  about  -2°C. 

5.  Both  the  initial  tangent  modulus  and  the  50% 
strength  modulus  are  found  to  be  independent  of 
strain  rate,  but  significantly  increase  with  the  de¬ 
crease  in  temperature  and  tend  to  decrease  with 
the  increase  in  density.  The  initial  tangent  modu¬ 
lus  found  in  this  study  is  about  17%  of  the 
Young’s  modulus  determined  by  Kaplar  (1963)  at 
1  °C  using  a  dynamic  method  and  83%  at  -10°C. 
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Figure  1.  Compacting  specimens  in  the 
Lucite  gang  mold. 
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Figure  9.  The  Model  1015  Biomat  ion  transient  recorder  and  a  Type  203  os¬ 
cilloscope. 
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Figure  10.  Typical  set  of  stress-strain 
curves  for  various  strain  rates  at  -5  °C. 
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Figure  II.  Typical  set  of  stress-strain  curves  for 
various  densities  at  a  strain  rate  of  IO  '*/s  and  a  tem¬ 
perature  of  -5  °C. 
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Figure  12.  Stress-strain  curves  for  various  densities  at  a 
strain  rate  of  10  '/s  and  a  temperature  of  -5°C. 
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Figure  13.  Typical  Fairbanks  silt  specimens  after  tension  tests. 
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Figure  17.  Peak  tensile  strength  as  a  function  of  time  to 
failure  for  the  specimens  with  various  densities  at  ~5  °C. 
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Figure  18.  Peak  tensile  strength  as  a  function  of  time  to 
failure  for  the  medium-density  samples  at  -2  and  -5°C. 
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Figure  19.  Plot  of  log  1  /a  vs  tm  for  the  medium-density 
samples  at  -2  and  - 5  °C. 
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Figure  22.  Peak  tensile  strength  as  a  function 
of  dry  density  for  different  machine  speeds 
at  -5  °C. 


Figure  20.  Log-log  plot  of  peak 
tensile  strength  vs  temperature  for 
a  strain  rate  of  IO  ‘/s. 


Figure  21.  Unfrozen  water  content  as  a  function 
of  temperature  for  Fairbanks  silt  with  a  water  con¬ 
tent  of  40.5%. 
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Figure  23.  Peak  strength  as  a  function  of 
volume  ice  content  for  different  machine 
speeds  at  -5  °C. 
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Figure  24.  Failure  strain  as  a  function  of  average  strain  rate  for  various  dry  densities  at  -5  °C. 


Figure  25.  The  variation  of  failure 
strain  with  temperature. 


Figure  26.  Comparison  between  the  tensile  failure  strain  and  the  compressive 
failure  strain. 


Measured  Strain, %  (neck  section ) 


Figure  33.  Comparison  between  the  measured  and 
the  corrected  strain  in  the  neck  section  of  speci¬ 
mens. 


Figure  34.  Comparison  between  the  a-t  curves  of 
the  coated  and  noncoated  specimens. 


Specimen 

no. 

W  ater 

content 

w 

(%) 

j 

density 

(>3, 
(8/cm  ) 

Volume 
ice  content 
i 

V 

Saturation 

degree 

S 

(%> 

Void 

ratio 

e 

Porosity 

n 

14-4 

38.8 

1.233 

0.454 

95.6 

1.174 

0.540 

15-1 

32.6 

1.367 

0.438 

98.4 

0.960 

0.4  90 

15-2 

31.5 

1.377 

0.424 

97.3 

0.946 

0.486 

15-3 

31. 1 

1.392 

0.423 

98.3 

0.925 

0.481 

15-4* 

30.4 

1.373 

0.407 

92.5 

0.952 

0.488 

16-1 

48.7 

1.108 

0.  549 

99.7 

1.419 

0.587 

16-2 

46.8 

1. 112 

0.528 

96.4 

1.410 

0.585 

16-3* 

47.6 

1. 120 

0.542 

99.9 

1.393 

0.582 

17-2 

34.6 

1.330 

0.455 

99.6 

1.015 

0.504 

17-4 

37.8 

1.268 

0.478 

98.5 

1. 114 

0.527 

18-1 

48.3 

1.112 

0.546 

99.6 

1.410 

0.585 

18-2 

51. 1 

1.098 

0.573 

100.0 

1.441 

0.590 

18-3 

48.9 

1.096 

0.546 

98.4 

1.445 

0.591 

19-1 

39.8 

1.243 

0.4  72 

99.5 

1.157 

0.536 

19-2 

40.9 

1.226 

0.481 

99.9 

1.185 

0.542 

19-3 

40.4 

1.243 

0.480 

100.0 

1.  155 

0.536 

19-4 

38.1 

1.274 

0.460 

99.8 

1.104 

0.525 

*  Specimen  coated  with  silicone. 
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Temperature 

0 

(°C) 

Unfrozen 
water  content 

w 

u 

(X) 

Tempe  rature 

6 

(°C) 

Unfrozen 
water  content 

w 

u 

(X) 

Temperature 

0 

(°C) 

Unfrozen 
water  content 

w 

u 

(X) 

-15.55 

3.06 

-18.45 

1.90 

-18.50 

2.17 

-14. 12 

3.18 

-15.44 

1.98 

-15.57 

2.21 

-10.85 

3.13 

-13.15 

2.06 

-13.39 

2.23 

-  8.95 

3.48 

-11.06 

2.13 

-11.09 

2.26 

-  6.86 

3.55 

-  8.40 

2.53 

-  8.48 

2.71 

-  4.96 

4.13 

-  7.36 

2.72 

-  7.44 

2.80 

-  3.86 

4.43 

-6.15 

2.80 

-  6.09 

3.03 

-  2.62 

4.90 

-  4.43 

3.19 

-  4.51 

3.56 

-  1.76 

5.55 

-  3.62 

3.65 

-  3.52 

3.86 

-  1.48 

6.42 

-  2.36 

4.31 

-  2.44 

4.45 

-  1.06 

7.83 

-  1.40 

6.25 

-  1.1 

5.57 

-  0.88 

7.95 

-  0.82 

6.78 

-  0.82 

6.90 

-  0.54 

8.12 

-  0.30 

13.77 

-  0.28 

13.85 

0.51 

0.38 

0.33 

0.28 


8.53 

9.64 

11.86 

16.31 


*  After  Tice  (unpublished) 
t  After  Xu  (unpublished) 
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A  facsimile  catalog  card  in  Library  of  Congress  MARC  format  is  repro¬ 
duced  below. 
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